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TECHNICAL MEMORANDUM 


PERFORMANCE ASSESSMENT OF LOW PRESSURE 
NUCLEAR THERMAL PROPULSION 


INTRODUCTION 


An increase in Isp for nuclear thermal propulsion systems is desirable for reducing the pro- 
pellant requirements and cost of future applications, such as the Mars transfer vehicle. Several 
previous design studies have suggested that the Isp could be increased substantially with hydrogen 
dissociation/recombination. Hydrogen molecules (H 2 ), at high temperatures and low pressures, will 
dissociate to monatomic hydrogen (H) (fig. 1). The reverse process (i.e., formation of H 2 from H) is 
exothermic. The exothermic energy increases the kinetic energy and, therefore, increases the Isp. 



TEMPERATURE (K) 


Figure 1. Mole fraction of hydrogen dissociated to monatomic hydrogen. 1 

The LPNTP system is expected to maximize the hydrogen dissociation/recombination and Isp 
by operating at high chamber temperatures and low chamber pressures. The process involves hydro- 
gen flow through a high-temperature, low-pressure fission reactor and out a nozzle. The high tem- 
perature (-3,000 K) of the hydrogen in the reactor is limited by the temperature limits of the reactor 
material. The chamber pressure is about 1 atm because lower pressures decrease the engines 
thrust-to-weight ratio below acceptable limits. This study assumes that hydrogen leaves the reactor 
and enters the nozzle at the 3,000 K equilibrium dissociation level. The computer model used begins 
with equilibrium at the nozzle inlet. 




Hydrogen dissociation in the reactor does not affect LPNTP performance like dissociation in 
traditional chemical propulsion systems, because energy from the reactor resupplies energy lost due to 
hydrogen dissociation. Recombination takes place in die nozzle due primarily to a drop in temperature as 
the Mach number increases. However, as the Mach number increases beyond the nozzle throat, the static 
pressure and density of the flow decreases which minimizes the recombination. The ideal LPNTP Isp at 
3,000 K and 10 psia is 1,160 s due to the added energy from fast recombination rates. The actual Isp 
depends on the finite kinetic reaction rates which affect the amount of monatomic hydrogen recombination 
before the flow exits the nozzle. 

An LPNTP system has other technical issues (e.g., flow instability and two-phase flow) besides 
hydrogen dissociation/recombination which affect the systems practicality. In this study, only the effects 
of hydrogen dissociation/recombination are examined. 


HYDROGEN DISSOCIATION/RECOMBINATION SENSITIVITIES 


A computer code was configured to model an LPNTP engine. The baseline nozzle design is shown 
in figure 2. The following model parameters were chosen based on a conventional conical nozzle design to 
keep the calculations simple for an initial study: 

• 5.2336-in throat radius 

• 27.6:1 nozzle contraction ratio 

• 200: 1 nozzle expansion ratio 

• Conical diverging nozzle section 

• Hydrogen propellant 

• 10 psia chamber pressure 

• 3,000 K chamber temperature 

• Nominal reaction rate coefficient 
Equations from Cohen & Westberg(C&W) 2 
k H = 3.2*10 15 

k m = (1.0*10 17 )*T-°- 6 

• All nozzle flow in one zone 

• Neglect boundary layer calculations and heat loss. 

Unless otherwise stated, all parameters are maintained at the baseline value for each computer run. 
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r =13.28 cm (5.23 in) 


Area* = 554.38cm 2 (85.93 in 2 ) 



Area = 15301.1 cm 2 
(2371.67 in 2 ) 
Area Ratio = 27.6 :1 


187.86 cm (73.96 in) 



Area = 110877 cm 
(17186 in 2 ) 

Area Ratio = 200 :1 


Figure 2. Baseline LPNTP nozzle design. 


KINETICS MODEL 


The two-dimensional kinetics (TDK) nozzle performance computer program 3 was used to 
determine the effect of various parameters on hydrogen dissociation/recombination and Isp. TDK 
simulates an invisid flow from the start of nozzle contraction to the nozzle exit. The one-dimensional 
equilibrium (ODE) option of TDK was used to calculate the upper Isp limit due to shifting molar species 
concentrations in the expansion process. The one-dimensional frozen (ODF) option of TDK was used to 
calculate the lower Isp limit due to frozen molar species concentrations in the expansion process. The one- 
dimensional kinetics (ODK) option of TDK accounts for the effects of finite chemical kinetics. Both ODE 
and ODK calculations provide data in the axial direction only, and the properties for each cross section are 
constant. The TDK option calculates changes in the radial direction and uses two third-body reactions for 
hydrogen (H 2 , H): 


H+H+H o H 2 +H 

(1) 

h+h+h 2 « h 2 +h 2 

(2) 


The reaction rate coefficient equations (£ Hz and kn) for each of the above reactions are in the 
Arrhenius form, which is written as: 


- 1,0008 

K = AT~ n e W . (3) 

Empirical data from C&W was used to determine the constants in the Arrhenius equation for both 
hydrogen third-body reactions. There is an uncertainty in k}\ and & H2- which affects the accuracy of the 

computed Isp. The uncertainty is due to the substantial variation in the results of hydrogen reaction rate 
experiments. Table 1 shows the forward kinetic and k^ equations obtained from C&W, from the 

national aerospace plane (NASP) rate constant committee of the NASP high-speed propulsion technology 
team, 4 and from the default equations used in the TDK. The C&W equations were used in this study. 
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Table 1. Kinetic reaction rate coefficient equations. 
H+H+M <=> H 2 +M 


Standard TDK 3 

M = At 

* = 6.4*10**(17)r**(-l) 


M = H 

k = 25**(Ar) 


M = H2 

* = 4**(Ar) 


Cohen and Westberc 2 

M = H 

*= 1 * 10 **( 15 ) 

Low 


* = 3.2*10**(15) 

Nominal 


*= 1.0* 10* *(16) 

High 

M = H2 

k = 5*10**(16)*7 1l, *(-0.6) 

Low 


* = 1 .0* 10**(17)*T**(-0.6) 

Nominal 


* = 2.0* 10**07)* r**(-0.6) 

High 

NASP Eauations 4 

M = H 

*= 1.5* 10**09)* r**(-1.0) 


m = h 2 

*= 1.8*10**(18)* P*(-1.0) 



Units = crn**6/(Mole**2*s) 


RESULTS 
Pressure Sensitivity 

This section shows the effects of chamber pressure ( Po ) on the Isp. The ODK and TDK Isp 
curves are close to equilibrium (ODE) at high Po, but approach frozen flow (ODF) at low Po (fig. 3). 
As Po decreases from 1,000 to 100 psia, both ODK Isp and TDK Isp increase. At 100 psia, figure 4 
shows the ODK Isp to be approximately 1,007 s, which is 27 s greater than the ODF Isp. From 100 
psia down to 10 psia, the ODK Isp and TDK Isp change is minimal. Below 10 psia, the ODK Isp 
increases like frozen flow due to a lower average molecular weight caused by greater dissociation 
and less recombination. Because Po below 1 atm was considered impractical, 10 psia is the lower 
limit of Po. At 10 psia, the ODK Isp is 1,005 s, but is only 14 s greater than the ODF Isp. 

The mole fraction of monatomic hydrogen (X H ) entering the nozzle is 17.3 percent with Po at 
10 psia and 5.9 percent with Po at 100 psia (fig. 5). The X H decreases in the nozzle due to mon- 
atomic hydrogen recombination. Inside the nozzle, there is 3-percent recombination at 100 psia and 
1.25-percent recombination at 10 psia, as shown in figure 5. Greater recombination from higher pres- 
sures is the reason why the difference between ODK Isp and ODF Isp is greater at 100 psia than at 
10 psia, as shown in figure 4. 
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Chamber Pressure (psia) 

Figure 3. Isp versus Po ( To = 3,000 K). 
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Figure 4. Isp versus Po (To = 3,000 K). 
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Figure 5. Xh versus axial nozzle position. 






The amount of one-dimensional Isp increase above frozen flow Isp (at 1,000 psia and 3,000 
K) due to hydrogen dissociation (lower molecular weights) and monatomic hydrogen recombination 
(added kinetic energy) is shown in figure 6, at various pressures. 



Figure 6. One-dimensional Isp increase from hydrogen dissociation and recombination 

versus chamber pressure. 


Temperature Sensitivity 

This section shows the effect on Isp by increasing the chamber temperature (To) to 3,200 K. 
Plots illustrating Isp versus Po at 3,200 K are shown in figures 7 and 8. The Isp is almost constant 
from 10 to 100 psia. At 100 psia, the ODK Isp is 1,066 s, which is 45 s greater than the ODF Isp. 
Thus, increasing To from 3,000 K to 3,200 K will increase the ODK Isp by approximately 60 s. 


1,500 
1,400 
_ 1,300 

f 1,200 
- 1,100 
1,000 
900 

1 10 100 1,000 10,000 

Chamber Pressure (psia) 

Figure 7. Isp versus Po (To = 3,200 K). 
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Figure 8. Isp versus Po ( To = 3,200 K). 


The effect of chamber temperature on the X H is shown in figure 5. Increasing To from 3,000 K 
to 3,200 K, at 10 psia, increases X H at the nozzle inlet to 28.6 percent. Results show 2-percent re- 
combination inside the nozzle. In addition, the exiting X H at 10 psia is higher at 3,200 K than at 
3,000 K. Greater X H , greater recombination, and higher To will all increase the Isp due to smaller 
average molecular weights and greater kinetic energy in the nozzle exhaust flow. 


Nozzle Contour Sensitivity 

The ODK results show X H at the nozzle exit to be 16.17 percent for the baseline nozzle con- 
tour. This section shows the sensitivity of the nozzle contour (within TDK limits) on monatomic 
hydrogen recombination. 

Table 2 displays the effects of nozzle geometry variations on X H . Results show that by vary- 
ing either the contraction ratio or the expansion ratio, there is no effect on the recombination. Holding 
the area ratios constant and lowering the nozzle converging angle from 30° to 10°, lowers the X^ at 
the nozzle exit to 15.9 percent, but increases the TDK Isp by only 3 s over the baseline. Lowering 
the nozzle diverging angle from 15° to 5° lowers the X H at the nozzle exit to 15.7 percent, and 
increases the TDK Isp 19 s above the baseline value. The higher Isp for small diverging cone angles 
is due to a greater axial component of exhaust velocity created by less divergence of the nozzle flow 
at the exit. Small converging and diverging cone angles require longer nozzles which may be 
impractical. For example, decreasing the diverging angle from 15° to 5° increases the total nozzle 
length from 300 in to over 800 in. 

Variation of the nozzle throat geometries is also considered because figure 5 indicates that 
most recombination occurs around the nozzle throat. Larger throat radii increase the recombination 
and TDK Isp (table 2). By increasing the throat radius from 5.23 to 10.0 in, the X H at the nozzle exit 
decreases to 15.6 percent, and the TDK Isp increases 6 s greater than the baseline TDK Isp. How- 
ever, increasing the throat radius from 5.23 to 10 in almost doubles the nozzle length, and increases 
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the thrust and mass flow rate by 365 percent. Other results show that the increased nozzle length increases 
propellant stay time in the nozzle, and indicate minimum changes in pressure and temperature at the nozzle 
throat compared with baseline LPNTP. Thus, some additional recombination could be caused by a longer 
stay time around the nozzle throat. 

The normalized wall radius of curvature downstream of the throat (RWTD), and the normalized 
wall radius of curvature upstream of the throat (RWTU) are displayed in figure 9. 


Table 2. Effects of nozzle geometry variations on recombination. 




Isp (s) 

Parameters 

Inlet (ODE) 

1.5 (ODE) 




ODK 

TDK 

Contraction Ratio 
low = 10 

17.30 



HI 



988 

nominal = 27.6 


N/A 




N/A 

988 

high = 50 

17.30 




WnSMm 


988 

Expansion Ratio 
low = 75 





16.17 


976 

nominal = 200 


N/A 


N/A 

16.17 

N/A 

988 

high = 500 





16.17 


995 


■ 







low =10° 
nominal = 30° 
hieh = 50° 

17.30 

17.30 

17.30 

N/A 




N/A 

991 

988 

988 

Diverging Angle 








low = 5° 



16.42 


15.74 


1,007 

17.30 

N/A 

16.41 


16.17 


988 

nominal = 15° 
high = 45° 

17.30 

16.41 


844 


17.30 


16.81 


■ 


982 

nominal = 5.23 

17.30 

N/A 

16.41 


■ 


988 

high = 10.0 

17.30 


16.05 




994 

RWTU 
low — 0.5 


16.64 

16.56 

16.33 


1,004 


nominal = 1.5 

N/A 

16.64 

16.41 

16.25 

N/A 


N/A 

high = 3.0 


16.64 

16.22 

16.01 


1,007 




16.64 

16.43 

16.28 


1,004 


nominal = 0.4 

N/A 

16.64 

16.41 

16.25 

N/A 

1,005 

N/A 

high = 1.5 


16.64 

16.40 

16.18 


1,005 


Contour 

Cone 

N/A 

16.64 

16.41 

16.25 

16.17 

1,005 

988 

Bell (100% length) 

16.64 

16.41 

16.29 

16.24 

1,004 

993 


Note: All cases were run at Po = 10 psia and To - 3,000 K. 


By doubling the RWTU, the X H at the 1.5 area ratio in the diverging section decreases from 16.25 
to 16.01 percent, and the ODK Isp increases from 1,005 to 1007 s. In addition, by increasing RWTD 
from 0.4 to 1.5, X H at the 1.5 area ratio in the diverging section of the nozzle decreases from 16.25 to 
16.18 percent, but the ODK Isp is almost constant. The effect due to changes in RWTU are greater than 
changes in RWTD because a larger RWTU increases both the converging nozzle length and the propellant 
stay time in the part of the nozzle which has higher pressures and flow densities. 
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<L 

Figure 9. Wall throat radius of curvature. 

The one-dimensional variations of pressure, temperature, and X H throughout the LPNTP 
nozzle are displayed in figures 10 and 11. Both figures show most recombination occurring around 
the throat. The change in X^ throughout the nozzle is very similar to changes in pressure. At the 
point were recombination stops, the pressure is close to 0 psia. However, the temperature continues 
to drop with no changes in X^. These results indicate that low static temperatures allow recombina- 
tion to begin, but low static pressures and flow densities limit the amount of recombination. 

The one-dimensional variations of pressure, temperature, and X^j throughout the baseline 
nozzle at 100-psia chamber pressure are shown in figures 12 and 13. As before, the recombination 
stops when the static pressure is close to 0 psia. There is greater recombination at higher pressures 
even though X H at the nozzle inlet is less. Thus, recombination appears greater in nozzle regions 
with high static pressures. The longer the stay time in these regions, the greater the amount of 
recombination. 
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-100 0 100 200 300 
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To = 5,400 °R 
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Figure 10. X H and P versus X ( Po = 10 psia). 
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The effect of a bell nozzle’s diverging contour on recombination was also examined. Table 2 
shows the effects of a bell nozzle contour on the X H and Isp. The bell nozzle has a slightly larger X^ 
at the nozzle exit, and an ODK Isp 1 s lower than the baseline conical nozzle. These results could be 
due to the bell nozzle having a greater diverging angle beyond the throat, which reduces the stay 
time in pressure regions needed for recombination. The bell nozzle has a greater TDK Isp because of 
less divergence losses at the nozzle exit. 


Kinetics Sensitivity 

The effects of hydrogen kinetic reaction rate uncertainty and key parameters in the Arrhenius 
equation were also analyzed. Although differences in reaction rate coefficients should not affect the 
results of a comparative study such as this, it will be an important aspect in ensuring TDK accuracy 
and Isp predictions. Previously shown in table 1 were three different sets of and % 2 equations. 

The effect of each equation on ODK Isp is shown in table 3. 

Table 3. Effects of various k H and £ H2 equations on ODK Isp. 

Hydrogen Reaction Rate ODK Isp (s) 

Coefficient Equations To = 3.000 K Tq = 3.2Q0-K 

C&W 1,005 1,068 

NASP 1,006 1,073 

Default TDK 1,008 1,076 

Results show a small difference in ODK Isp between the three sets of £ H and equations. 

Each set of equations was determined from experimental reaction rates. The uncertainty for each set 
of ifc H and k^ 2 equations is due to either scatter in experimental data or experiment type. 

Table 1 shows the C&W and *h 2 equations used to determine the impact of uncertainty on 
Isp. The nominal equations were recommended by C&W. The high and low equations of k H and &h 2 

at 3,000 K represent the fastest and slowest rate limits, respectively, based on the C&W level of 
uncertainty. The effects on Isp caused by C&W uncertainty are shown in table 4. There is an Isp in- 
crease when either/or both Jt H and Jfc H2 become larger. Note, there is a slightly greater effect on Isp 
when varying k H instead of k Hj due to the greater uncertainty of the monatomic equation. The ODK 

Isp goes from 999 s with the slowest rates to 1,014 s with the fastest rates. Table 5 shows similar 
results, but with To equal to 3,200 K. At 3,200 K, the ODK Isp goes from 1,056 s with the slowest 
rates to 1,085 s with the fastest rates. Higher temperatures decrease the Isp accuracy because of 
greater uncertainty. The maximum variation in Isp due to uncertainty is significant and should not be 
ignored. 

The A and N parameters in the Arrhenius equation were examined to determine the impact on 
the Isp. By multiplying the nominal value of A in and k^ 2 by various assumed multipliers, the ODK 

Isp will increase as shown in figure 14. Rates which are 10 times faster than nominal will increase 
the ODK Isp by approximately 30 s. The TDK Isp is more sensitive to variations of N in the 


11 



Arrhenius equation for than k H , as shown in figures 15 and 16. Decreasing A for k U2 from 1.0 to 

0.2, increases the TDK Isp by about 44 s. Isp effects when varying both A and N for each computer 
run were not examined. 

Table 4. Isp variations with reaction rate accuracy (To = 3,000 K). 

H+H+M <=> H 2 +M 



M = H 

M = H 2 

Low 

Nominal 

High 

Low 

999* 

1,002 

1,009 


(982)** 

(985) 

(992) 

Nominal 

1,001 

1,005 

1,011 


(985) 

(988) 

(994) 

High 

1,007 

1,009 

1,014 


(990) 

(993) 

(997) 


Isp (ODE) = 1,159 s *ODK Isp (s) 

Isp (ODF) = 991 s **TDK Isp (s) 

To = 3,000 K 
Po =10 psia 
Ae/A* = 200 

Table 5. Isp variations with reaction rate accuracy (To - 3,200 K). 

H+H+M <=>H 2 +M 


M = H 2 

M = H 

Low 

Nominal 

High 

Low 

1,056* 

(1,039)** 

1,064 

(1,047) 

1,079 

(1,062) 

Nominal 

1,061 

(1,043) 

1,068 

(1,050) 

1,081 

(1,064) 

High 

1,068 

(1,051) 

1,074 

(1,056) 

1,085 

(1,067) 


Isp (ODE) = 1,297 s 
Isp (ODF) = 1,041 s 
To = 3,200 K 
Po = 10 psia 
Ae/A * = 200 


*ODK Isp (s) 
**TDK Isp (s) 
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Kinetic Reaction Rate Multiplier 


Figure 14. Isp dependence on kinetic reaction rate magnitude. 



Figure 15. Isp dependence on the Arrhenius temperature exponent for H. 



Temperature Exponent in the Arrhenius Equation 

Figure 16. Isp dependence on the Arrhenius temperature exponent for H 2 _ 
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Best Potential LPNTP Design 


Based on the parameters examined, the optimum value for each parameter and the fastest 
reaction rates from C&W were input into the TDK model. The following input parameters are dif- 
ferent from the baseline: 

• Converging angle = 5° 

• Contraction ratio = 7.41:1 

• RWTU = 3.0 

• RWTD = 1.5 

• Bell nozzle with diverging section 
78-percent baseline diverging nozzle length 

• p 0 = 44 psia 

• To = 3,200 K 

• * H = 1.0*10 16 

• *h 2 = (2.0*10 17 )*T-°- 6 . 

The total nozzle length is the same as the baseline nozzle length. The ODK Isp result is 
1,080 s. In addition, the ODK thrust is 6,845 lbf and the mass flow rate is 6.338 lbm/s. The ODK Isp 
result is slightly less than the largest ODK Isp from table 5, which uses the fastest rates with base- 
line nozzle geometry. This result shows minimum impact on ODK Isp from optimum nozzle geometry. 


CONCLUSIONS 


TDK Isp is close to equilibrium flow at high Po and close to frozen flow at low Po. With To at 
3,000 K, the maximum TDK Isp and maximum difference between ODK Isp and ODF Isp occurs at 
Po -100 psia, due to increased monatomic hydrogen recombination. The Isp is -15 s higher at a Po of 
100 psia than at 1,000 psia. The optimum Po indicates a greater impact on Isp from recombination 
than from lower average molecular weights. A Po of 3,200 K increases both the dissociation 
(decrease in the average hydrogen molecular weight) and the recombination of the flow. These 
characteristics boost the Isp above the value that would be obtained by only considering temperature 
effects. 


Examination of various nozzle geometries showed a minimum impact on recombination and 
Isp. The significant changes that did occur were related to a more efficient use of two-dimensional 
nozzle flow. Most recombination occurs around the nozzle throat due to low static temperatures 
which foster recombination. Low static pressures and flow densities limit the amount of recombina- 
tion downstream of the throat. 

Finally, the uncertainty of the and & H2 , at high To, has a significant impact on the Isp pre- 
dictions. There is a large difference between ODE and ODF Isp’s at a lower Po, and the ODK Isp 
could be anywhere between, depending on the assumed reaction rate. The slow and fc H2 used in 

this study cause the ODK Isp to approach the ODF Isp at low static pressures. 
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RECOMMENDATIONS 


Since the LPNTP Isp predictions were closer to frozen flow Isp than ideal equilibrium Isp at 
low Po, more extensive calculations, using a TDK model that accounts for boundary layer effects, 
and nonadiabatic flow are not recommended at this time. 

Faster kinetic reaction rates will increase both the monatomic hydrogen recombination and 
the Isp. Because there is an uncertainty with the high temperature values of and k^ 2 , accurate 
hydrogen dissociation/recombination bench tests are recommended to determine if the actual reaction 
rates are different than the published reaction rates used in this study. High temperature subscaled 
nozzle tests with load cells and flowmeters are recommended to determine actual Isp, because 
actual nozzle flow characteristics (e.g., boundary layer effects) might affect the amount of hydrogen 
dissociation/recombination. 

Based on the results of this study, the chamber pressure recommended is 190 psia for a 3,000 
K chamber temperature. This chamber pressure is above hydrogen’s critical pressure (188 psia) to 
minimize two-phase flow problems in the feed system (e.g., flow oscillations). Operating at 190 psia 
has an Isp slightly less than at 100 psia, but produces a higher thrust than at 10 psia and has a 
higher Isp than at 1,000 psia. In addition, the nozzle should have a bell diverging section to shorten 
the nozzle length and minimize divergence losses. The nozzle’s converging/throat section should be 
designed to maximize the recombination. 

Finally, TDK only uses conical converging nozzle sections. Other kinetic computer models 
which can vary the converging nozzle contour should be investigated to determine the effects on Isp. 
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